Equation of state and phonon frequency calculations of diamond at high pressures. 
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The pressure-volume relationship and the zone-center optical phonon frequency of cubic diamond 
at pressures up to 600 GPa have been calculated based on Density Functional Theory within the 
Local Density Approximation and the Generalized Gradient Approximation. Three different ap- 
proaches, viz. a pseudopotential method applied in the basis of plane waves, an all-electron method 
relying on Augmented Plane Waves plus Local Orbitals, and an intermediate approach implemented 
in the basis of Projector Augmented Waves have been used. All these methods and approximations 
yield consistent results for the pressure derivative of the bulk modulus and the volume dependence 
of the mode Griineisen parameter of diamond. The results are at variance with recent precise mea- 
surements up to 140 GPa. Possible implications for the experimental pressure determination based 
on the ruby luminescence method are discussed. 

PACS numbers: PACS: 71.15.Nc, 63.20.-e, 62.50,+p 



I. INTRODUCTION 

Diamond is the archetype of the covalently bonded, 
tetrahedrally coordinated insulators. It's extreme hard- 
ness is highly valued in technology and is also exploited in 
high pressure research when using the diamond anvil cell. 
The elastic properties of diamond near ambient condi- 
tions are well characterized through ultrasonic and Bril- 
louin techniques. However, as a consequence of the small 
compressibility, the changes in elastic properties under 
high hydrostatic pressure are not well confined experi- 
mentally. This, for instance, applies to the variation of 
the bulk modulus B with pressure P, B' = (dB /dP)p—o, 
a basic parameter in the equation-of-state (EOS) mod- 
elling. A property closely related to the compression 
behavior of diamond is the pressure-induced frequency 
shift of the threefold degenerate F2 S zone-center optical 
phonon mode; its shift with pressure provides an approx- 
imate measure of the change in relative density, because 
the mode Griineisen parameter is close to onei 1 * 2 * 3 *^ 

High-pressure x-ray diffraction experiments** 6 * 7 - yield 
the ambient-pressure bulk modulus Bq in good agree- 
ment with acoustic measurements ^a2ii**iAii2 The analysis 
of diffraction data is usually based on adopting B' = 4.0 
obtained from ultrasonic measurements up to 0.2 GPa. 8 
The only exception is the recent diffraction study of dia- 
mond to 140 GPa by Occelli et al.^ who report B' = 3.0, 
at variance with the ultrasonic measurement. It was sub- 
sequently argued 13 that the ruby pressure calibration** 
employed in the diffraction work of Ref. @ may need 
a revision. If so, the pressure shift of the F2 9 phonon 
frequency, also studied by Occelli et al. up to 140 GPa, 
should be affected in a similar manner. The F2 g phonon 
mode behavior at high pressures was frequently stud- 
ied by Raman spectroscopy 1 ! 2 ! 3 ! 4 ! 5 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 and its 
possible role in pressure calibration was addressed early 
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We report the calculation of the EOS and optical 
phonon frequency of diamond at high pressures within 
Density Functional Theory (DFT). Extensive theoretical 



work on diamond under pressure by ah initio methods 
has addressed changes in bonding, elasticity, lattice dy- 
namics, thermodynamical properties, phase stability, and 
electronic excitations, see , e.g., Refs. 322ll23ll24ll25ll2"ri 
W^^^^^^^^^^k?,1i?^?^^^^ . Here, 
we are interested in a specific question: What are, 
based on different implementations and approximations 
of DFT, the constraints on B' of diamond and on the 
nonlinear pressure shift of the optical phonon frequency. 
The calculated values of B' reported in the literature 
scatter by about 25%, spanning a range similar to that 
of the experimental results. Calculated pressure effects 
on the optical phonon frequency were reported in Refs. 
Q and |23 • Results of a more recent calculation 41 - show 
some disagreement with Refs. Q and j^. In view of 
new experimental results 7 - and the - at least apparent - 
uncertainties in the previous theoretical predictions we 
considered it worthwhile to revisit the calculation of the 
EOS and optical phonon frequency of diamond under 
pressure, combined with accurate procedures to extract 
the parameters of interest. 



II. DETAILS OF THE CALCULATIONS 

The total energy calculations performed in this work 
are based on DFT 43 within, on one hand, the Local Den- 
sity Approximation 4 *^ (LDA) and, on the other hand, the 
Generalized Gradient Approximation*^*** 6 - (GGA). In or- 
der to grasp the uncertainties consequent to the choice 
of the computational method and of its inherent assump- 
tions, we are using simultaneously three different ap- 
proaches, viz. the pseudopotential method applied in 
the basis of plane waves*^** 4 ***^ 9 * 5 *^ (PW), an all-electron 
method relying on the Augmented Plane Waves plus Lo- 
cal Orbitalo 51 ! 52 ! 53 (APW+lo), and an intermediate ap- 
proach implemented in the basis of Projector Augmented 
Wave o 54 ! 55 (PAW); the latter approach*! 6 - treats the va- 
lence states as part of an all-electron problem and de- 
scribes them by all-electron wave- functions. 
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The pseudopotential employed is the "dual-space" sep- 
arable pseudopotential of Hartwigsen, Goedecker, and 
Hutted and its GGA counterpart constructed by X. 
Gonz o 47 i 58 ; the PAW potentials were constructed by G. 
Kresse and denoted as C_h in Ref. 54]. All methods 
take into account scalar-relativistic corrections (though 
no meaningful contribution expected in diamond), either 
explicitly (APW+lo) or through the construction of the 
(pseudo-)potentials. 

The numerical convergence of all three methods, with 
respect to the size of the basis set and k-space sampling 
was thoroughly tested. The plane-wave cutoffs of 75 
Hartree and 1500 eV were applied with, respectively, the 
PW and the PAW-basis. These rather large values are re- 
quired for getting reliable results for the pressure which 
is calculated analytically, in both approaches, with the 
aid of the stress theoremiSiSi In the APW+lo method a 
muffin-tin radius of Rmt = 1.2 a.u. was used. The plane- 
wave expansion was limited by Rmt x K max = 9 and the 
charge density Fourier expansion by G max — 20 Ry 1 ^ 2 . 
In the APW+lo method the P(V) results were generated 
from E tot (V) fitted as described below. 

The integration over the Brillouin zone is performed 
using the same set of 28 k-points in the irreducible wedge 
in both plane-wave based approaches; the set was gener- 
ated by the "special points" approach- 1 using a 6 x 6 x 6 
mesh with 4 different fractional shifts. In the APW+lo 
approach the tetrahedron integration was employed^* 2 ^ 
based on a uniform, T-centered 10 x 10 x 10 mesh of 47 
k-points. 

The frozen-phonon approach, which consists in the 
evaluation of the total energy E tot of the crystal with 
frozen-in atomic displacements, was applied (in the two 
plane-wave based methods only) with the same cutoffs 
and with the same k-point mesh (6x6x6 + 4 shifts) 
which, due to the lowered symmetry, results in a 91 k- 
points set. Small displacements u(l) = ±(u,u,u) and 
u(2) — +(u,m, u) with u/a — 0.002 were applied to 
the two atoms of the basis - thus either compressing or 
stretching the C(l) - C(2) bond, like in the zone-center 
F2 g mode - and the two values obtained for Ai? tot [viz. 
the AE tot (outward) and Ai? tot (inward)] were averaged. 
The eigenfrequency is then found from the expression for 
the energy of a harmonic oscillator 



AE tot = l -Mu? [|w(l)| 2 + |u(2)| 2 



(1) 



where M is the atomic mass. The above-mentioned 
averaging procedure of the AE tot (outward) and 
AE tot (inward) eliminates the cubic contribution to 
AE tot . It turns out that the remaining quartic anhar- 
monicity is small at u/a — 0.002 and contributes to the 
uncertainty of the resulting eigenfrcquencies by less than 
0.1 cm -1 (checked by repeating the same calculations 
with displacements w/a=0.0015 and 0.003). 

Total energies, pressures, and phonon frequencies 
were calculated at 12 different volumes ranging from 
6.35 A 3 /atom down to 3.35 A 3 /atom, i.e., for about 10% 
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FIG. 1: Pressure- volume results for diamond from plane wave 
(PW, triangles up), projector augmented wave (PAW, trian- 
gles down), and all-electron (APW+lo, diamonds) calcula- 
tions. Solid and dashed lines represent fits of Eq. |5] to the 
combined GGA (filled symbols) and LDA (open symbols) re- 
sults, respectively. The inset offers a zoomed view of pressure 
differences relative to the average P(V) relation for GGA cal- 
culations, including the differences for experimental data of 
Ref. p] and their extrapolation. The volume dependence of 
the zero-point pressure P ZP (V) is also shown in the inset. 



volume expansion to 40% compression relative to the 
experimental equilibrium volume of 5.6725 A 3 /atom at 
300 K. — The lower volume limit corresponds to roughly 
600 GPa maximum pressure. 



III. PRESSURE- VOLUME RELATIONSHIP 

A. Analytical form 

Figure ^ shows calculated pressures as a function of 
volume. An excellent consistency of the data points ob- 
tained by the three different methods is immediately ap- 
parent. Also, the present results substantiate the claims 
made in Ref. [55j concerning the equivalence of the 
PAW approach with all-electron methods such as that 
of Ref. IH3. 

The proper analysis of the P(V) results yields the 
equilibrium properties Vq, Bq = —Vo(dP/dV)Q, and 
B' = (dB/dP) = (dlnB/cilnV^o (throughout this 
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FIG. 2: Bulk modulus Bo of diamond obtained from dif- 
ferent calculations plotted against the calculated equilibrium 
volume Vo- Symbols for the present GGA and LDA results 
as m Fig. □ Other calculated values of B |24l2^2al3al35l| 
are indicated by open circles (the attached letters refer to the 
first-author name). The experimental values of Bo originate 
from ultrasonic measurements^ and Brillouin scattering^- 
The solid (dashed) line refers to the volume dependence of 
the bulk modulus for combined GGA (LDA) results. The 
inset shows the pressure derivative of the bulk modulus ver- 
sus atomic volume. The calculated B' values and the two 
experimental ones are represented by symbols. 



Here, Eq is the energy at V = Vq. Equation [5] was used 
to fit the P(V) results obtained via the stress theorem 
(PW and PAW methods), while Eq. served to deter- 
mine the EOS parameters from total energies (APW+lo 
method) . For the PW and PAW methods it was carefully 
checked that the directly calculated P(V) data are con- 
sistent with those obtained by differentiation of E(V). 
For all sets of calculated P(V) or E(V) data the rms 
deviations were less than 0.04 GPa or 0.05 meV, respec- 
tively. 

Aleksandrov et al~ pointed out that a quadratic de- 
pendence of pressure on change in relative density 
Ap/p = Vq/V — 1, written as 



P{P) = B { 



Ap 

Po 



B' -lAp 



Po 



(4) 



applies in the case of diamond for P < 100 GPa. This 
expression is found to fit the calculated results up to 
600 GPa quite well. It is only in the statistical sense 
that Eq.0]is slightly inferior to Eq.[3 i.e., the rms devia- 
tions are larger by a factor of two, but this is not relevant 
for the present discussion. The bulk modulus, obtained 
by differentiation of Eq.^Jwith respect to normalized vol- 
ume, is given by 



B(p) = B 



P_ 

Po 



1 + (B' - 1) 



Ap 

Po 



(5) 



Equation[S]corresponds to a quadratic dependence of the 
bulk modulus on relative change in density. 



manuscript the subscript 'zero' refers to zero (ambient) 
pressure or equilibrium volume, either calculated or ex- 
perimental depending on context). One also needs to 
identify an analytical form of the P(V) behavior which, 
using the above parameters, fits the calculated pressures 
over the full volume range. 

We have tried several of the common equation-of-state 
forms iSiSSiSLS The best analytical expression, in the 
least-squares sense, was identified as 

P(V) = WoX~ n (l - X) expfa (1 - X)} (2) 

where X = (V/V a ) 1/3 , rj = 3B /2 + l/2-n, and n = 7/2. 
With n = 7/2 Eq. H is a blend of the Vinel£L (or 
Rydberg)£2i2i form and the Holzapfel^ expression, for 
which n = 2 and n = 5, respectively. Actually, the pa- 
rameter n varies between 3.4 and 3.6 for individual sets 
of calculated P(V) data of diamond. Without loss of sig- 
nificant digits in the fitted parameters of interest, we fix 
n at a value of 7/2; in this case the energy versus volume 
relation, obtained by integration of Eq. EJ can be written 
without invoking special functions other than the error 
function: 

E(V) = E + 9B Vo[f(V)-f(V )]en>(v)/Vv (3) 
f(V) = ^(277 + l)erf(07A 2 ) + 2^exp(-r;X)/X 2 . 



B. Calculated EOS parameters 

The obtained pairs of (Vo,i?o) values are displayed in 
Fig. Within a given approximation for the exchange- 
correlation functional (GGA or LDA), the three methods 
employed here (PW, PAW, APW+lo) yield results for 
Vo and Bq which are in very good agreement with each 
other. The choice of the exchange-correlation functional, 
however, leads to differences of 4% in Vo and 5% in B , 
reflecting in part the well-known overbinding of the LDA. 
All B' Q values (cf. inset to Fig. [5J fall into a narrow range 
between 3.6 and 3.7 , the difference between averaged 
LDA and GGA results being close to 1%. 

Since the calculated EOS parameters clearly split into 
two groups, i.e. the GGA and LDA results, with nearly 
identical parameter values within each group, we have 
combined the P(V) points for GGA and LDA, respec- 
tively, to obtain the GGA and LDA parameters Vq, Bo, 
and B' Q listed in Table [I] The corresponding P(V) re- 
lations are shown by solid and broken lines in Fig. ^ 
where the inset illustrates deviations of individual calcu- 
lated results from the average curves. Similarly, the solid 
and broken lines in Fig. [21 and its inset show the corre- 
sponding volume dependences of the bulk modulus and 
its pressure derivative. 
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TABLE I: Calculated equation-of-state [Eq. |5j and optical phonon parameters of diamond. Selected experimental results are 
listed in the lower part of the Table. The static-lattice entries correspond to switching off (at T = K) the effect of zero 
point motion. The calculated phonon pressure coefficients are obtained via (dui/dP)o = 70W0/B0 with parameters given in the 
same row of the Table. Likewise, the experimental mode Griineisen parameters were obtained from the experimental pressure 
coefficients unless noted otherwise. In the row marked GGA&LDA both the ambient pressure volume and bulk modulus are 
normalized to one. 



Method 


Vo 


Bo 


B' 


11 


wo 


70 


7o 


(duj/dP) 




A 3 


GPa 






cm -1 






cm" 1 /GPa 


All LDA 


5.510(5) 


465(3) 


3.63(3) 


7/2 


1322(2) 


1.003(3) 


0.79(5) 


2.87 


All GGA 


5.697(4) 


433(2) 


3.67(3) 


7/2 


1290(2) 


0.995(3) 


0.80(5) 


3.00 


GGA&LDA normalized 


1 


1 


3.65(5) 


7/2 


1 


1.000(5) 


0.80(5) 




Expt. T=300 K 


5.6725° 


442" 


4.0(5)" 




1332.5 


0.962(15) c 




2.90(5)° 






444.8(8) e 






1332.40(5) e 


1.00(3/ 




3.00(10/ 




5.674(1)° 


446(1) 9 


3.0(1) 9 


2 


1333 9 


7=0.97 9 




2.83 9 


Expt. T^OK 


5.6707° 


445* 






1332.70(3) 








Static lattice 


5.5886 ft 


469* 
















5.6122 3 


462* 















°Ref. pi . 

"Ref. pj; ultrasonic experiments up to 0.2 GPa. 
^Experime ntal mode G riineisen parameters reported in the liter- 
ature Ill2l3l4l5ll5ll7ll8l vary between 0.90 to 1.0 6. 
"'Best ' va lue in the literature according to Ref. I'2(l . 
e Ref. Il ll : Brillouin and Raman scattering. 

^Value obtained for a revised linear pressure coefficient of the ruby 
line shift as explained in the text. 

9 Ref. ; x-ray diffraction ( Vinet fit) and Raman scattering up to 
140 GPa. The mode Griineisen parameter 7 is the average for the 
pressur e ra nge 0—140 GPa. 

h Ref. 1711 : the classical limit in path integral Monte Carlo simu- 
lations. 

'Estimated value, see text. 

3 Ref. ll ll : based on extrapolation of isotope effects. 



The larger Bo value for LDA compared to GGA cor- 
relates with the smaller equilibrium volume. However, 
the B(V) curve for LDA falls below that for GGA (see 
Fig- 13 1 • This also applies to Bo values obtained in other 
self-consistent LDA calculations for diamond2LS2i^iS 
(represented by open circles in Fig. |2J . 



mimic the zero-point effects. 

The zero-point vibrational pressure P zp , i.e., the iso- 
choric change in pressure when zero-point motion is 
switched on, is about 6 GPa at the static-lattice value 
of Vo^i An approximate relation for the volume depen- 
dence of P zp is 



C. Zero-point motion effects 

The calculated results for Vo and Bo should be com- 
pared to experimental properties after 'correction' for vi- 
brational effects. Thermal effects are almost negligible 
at 300 K because of the high Debye temperature of di- 
amond (~2000 K), but zero-point motion needs to be 
considered. According to Monte Carlo simulations^ the 
static-lattice equilibrium volume for diamond is about 
5.59 A 3 . All the data for V shown in Fig. fall into a 
range of ±2% around this value. The static-lattice bulk 
modulus at V — 5.59 A 3 is estimated to be 469 GPa 
(this follows from the experimental values of Vo and Bo 
and the calculated B' ). Thus, relative to static-lattice 
properties, the GGA results cannot be considered supe- 
rior to the LDA ones; they appear better when compared 
to experimental data because the GGA errors happen to 



dV 



where AV ZP is the isobaric volume expansion due to 
zero-point motion. For getting AV ZP (V) within the 
quasi-harmonic thermodynamics, we simply refer to the 
calculated pressure dependence of the zero-point expan- 
sion presented in Fig. 10 of Ref. The corresponding 
volume dependence of P zp is displayed in the inset of 
Fig. HI P zp increases by about 6 GPa at 600 GPa. 
Adding such small changes in P zp to the calculated 
P{V) results leads to an increase of the B values by 
about 1%, but would not affect the calculated B' Q values 
(within the estimated uncertainty of the optimum value 
given below). 
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D. The essence of the EOS calculations 

In Fig. |3 the calculated P(V) relations are plot- 
ted in reduced coordinates, i.e. pressure normalized by 
Bq (Ap/po) as a function of A p/po- In this representation 
small differences in ambient-pressure volumes and bulk 
moduli are suppressed and a linear slope corresponds to 
a quadratic dependence of pressure on change in relative 
density (Eq.^J). The LDA and GGA results hardly differ 
in slope. In this way, Fig.[3]illustrates the main informa- 
tion we extract from our EOS calculations for diamond: 

(1) The value of B' does not depend much on crystal 
potential and basis set issues and on approximations for 
the exchange-correlation functional. Our average value 
is 

B' a = 3.65 ±0.05, 

where the small uncertainty reflects the scatter for the 
different methods of calculation. 

(2) Within the volume range considered here, the most 
appropriate 3-parameter analytical form for the P{V) re- 
lation of diamond (Eq. or is transferable between 
LDA and GGA solutions; only Vq and Bo need to be 
adjusted. 

(3) Inserting the calculated B' value and the experi- 
mental data for V and B (see Table[j} into Eq.[U(Eq.|lJ| 
is considered to yield, on the basis of this work, the op- 
timum representation of the EOS of diamond at 300 K. 
Actually, at the experimental Vq the GGA bulk modulus 
(calculated value of 440 GPa plus 1% correction for zero- 
point motion effects) happens to be very close to the ex- 
perimental Bq, and the calculated pressure at the exper- 
imental Vq is only 2 GPa. Therefore, the optimum EOS 
corresponds to the P{V) relation obtained from GGA 
calculations, with only a small correction applied. 



1.4 



1.0 



I I I I I 
Diamond / 

1 ' 1 III \-J 1 1 \A / 

/ 


/ GGA 

X 


B'=4 / 

7 X 

/ / 






r . 

-'' 






// fx' 




/ / 

/ X / 1 


11 GPa 


/ X X ' 

7 X 140 GPa 




.X/'Experiment B'=3 








i i i i i 





0.0 0.1 0.2 0.3 0.4 

Relative density change Ap/p 



FIG. 3: Calculated pressure-volume relations of diamond 
in reduced coordinates and using Vo, Bo, and B' from Ta- 
bleQ Experimental data of Ref. [7J are shown for comparison. 
These experiments cover a pressure range of about 140 GPa. 
The line labelled B' — 4 represents the EOS with parameters 
Bo and B' from ultrasonic experiments 0. 

difference is illustrated in Fig. El At nominally 140 GPa, 
about the maximum pressure reached in the experiments, 
the difference between our calculated B' () value and that 
Ref. translates to a pressure difference of about +11 
GPa. This means that either repulsion is slightly overes- 
timated in the calculations (independent of the exchange- 
correlation functional) or that the experimental data suf- 
fer from systematic errors, a combination of these effects 
not being ruled out. We will return to this issue below. 



E. Comparison with other results for B' 



IV. OPTICAL PHONON FREQUENCY 



Calculated values of B' reported in the literature are 
3.54 IH, 3.6, 4.5 |3, 3.24 |3, 3.5 |H and 3.97 [Til- In 
some cases, the differences to our result are small. Larger 
deviations may in part result from the procedure used to 
extract the EOS parameters (e.g., fit of a Murnaghan 
equation which is inadequate). 

The experimental value B' =4 ± 0.5 (Ref. @) de- 
rived from sound speed measurements at low pressures 
(0.2 GPa) is not sufficiently accurate to test the calcu- 
lated result. 

The only other experimental value of B' stems from 
the recent x-ray diffraction experiments up to 140 GPa. 7 
The EOS data appear to be of high quality. They were 
measured using helium as a pressure medium which is 
considered to provide almost hydrostatic conditions in 
DAC experiments. The obtained value of B' Q = 3.0(1) 
is inconsistent with the ultrasonic result and also sig- 
nificantly smaller than the calculated value. The latter 



A. Calculated phonon frequency versus volume 

Figure 2] shows calculated F2 ff phonon frequencies as a 
function of volume. At a given volume, the frequencies 
of all calculations are quite consistent with each other. 
The frequencies are too low by about 3% if compared to 
the experimental phonon frequency at ambient pressure. 
This small discrepancy does not originate from anhar- 
monicity connected to the frozen-phonon displacement. 

The volume dependence of phonon frequencies is usu- 
ally characterized by the simple scaling law 



«>(V) 

LUQ 



(6) 



which assumes that the mode Griineisen parameter 7 = 
— dhiuj / d\nV is independent of volume. The frozen- 
phonon results indicate some volume dependence of 7, 
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FIG. 4: Calculated zone-center optical phonon frequency of 
diamond as a function of atomic volume. The solid line guides 
through the combined GGA results. A Griineisen relation 
Eq.|S|with parameters loq = 1333 cm -1 (experimental phonon 
frequency) and 70 = 1 is indicated for comparison (dashed 
line) . Differences of individual calculated results with respect 
to the 'combined GGA' curve are shown in an inset. The 
second inset illustrates the calculated volume dependence of 
the mode Griineisen parameter (for combined GGA and LDA 
results, respectively) . 



an observation which was also noted by NielseniSi The 
expression 



Auj(V) 



la 



"7o 



1 



(7) 



with Au)(V) — u>(V) — luq yields a slightly better match 
of our frozen-phonon results (3 cm -1 rms deviation or 
better for individual sets). Within standard deviations, 
the parameters 70 and 7q for the combined GGA and 
LDA results are identical (see Table QJ. The upper inset 
to Fig. 0| illustrates the small volume dependence of the 
mode Griineisen parameter resulting from "f' < 1. 

With p/po ~ LojujQ (7 w 1) it follows from Eq. [5] 
that the ratio of normalized bulk modulus to normal- 
ized phonon frequency is approximately linear in relative 
density and the slope is B' — 1: 



'(HY as 1 + (So - 1) — 
B \ujq J p 



(8) 



frequency could come handy as a simple relation between 
an elastic and a dynamical property of diamond. 



B. Calculated vs. 



experimental mode Griineisen 
parameter 



This scaling between bulk modulus and optical phonon 



To compare the calculated phonon results with ex- 
perimental data, we first consider the measured linear 
pressure coefficient of the phonon frequency. Accord- 
ing to Schiferl et aZ.jSi the 'best' experimental value in 
the literature, taken as the average of Refs. |3 ll5ll9| . 
should be {dtu/dP) = 2.90 ± 0.05 cm-yGPa ( 7o « 
0.96). This value originates from experiments which em- 
ployed the ruby pressure calibration with a linear co- 
efficient A = (dP/d\n\) = 1905(10) GPa for the Rl 
line wavelength (A) shift near zero pressure as deter- 
mined by Piermarini et al£m Recent high-precision mea- 
surements of the ruby line shift up to 1 GPa- 3 " and a 
rcintcrpretation 13 of Rl line shift data measured up to 
~20 GPa— i24 indicate that the parameter A is smaller 
(A = 1820 ± 30 GPa) compared to the previously ac- 
cepted value. This revision leads to a corresponding in- 
crease of the phonon pressure coefficient. In this context 
it is helpful that the ratio of the diamond phonon fre- 
quency shift to ruby wavelength shift is explicitly given 
in Refs. and 0- The results, 7.94(10) cm" 1 /™ and 
7.85(6) cm -1 /™, agree quite well with each other. Tak- 
ing the average value (7.9(1) cm _1 /nm) in combination 
with A = 1820(30) GPa, the corrected value from DAC 
experiments should be (du>/dP)o = 3.0(1) cm /GPa. 
Within experimental uncertainties, this linear pressure 
coefficient agrees with results obtained by methods which 
do not involve any ruby calibrationjii^ in particular that 
of Whalley et al. 2 With B =444.8 GPa the correspond- 
ing mode Griineisen parameter becomes 70 = 1.00(3), 
the error being mainly due to the uncertainty of the ruby 
coefficient A. We note the excellent agreement of the cor- 
rected experimental value for 70 with the results obtained 
within the LDA and GGA approximations (see Table [I]). 

Occelli et al. report an average mode Griineisen pa- 
rameter of 7 = 0.97 for the range 0-140 GPa. It should 
be pointed out that their value of 7 is independent of 
the pressure scale used in the experiments because the 
volume was measured directly. Taking into account the 
predicted volume dependence 7 (V) (see inset to Fig. @}, 
the calculations are fully consistent with the experimen- 
tal volume dependence of the phonon frequency. This is 
no longer the case if we turn to the pressure dependence. 



C. The pressure — phonon-frequency relationship 

One can solve Eq. [7| for the relative volume, insert 
the corresponding expression into Eq. [2] (Eq. 0J, and 
use the 'GGA' and 'LDA' parameters from Table [I] to 
obtain an analytical form for pressure as a function of 
phonon frequency. On the other hand, in the analysis 
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TABLE II: Calculated zone-center optical phonon properties 
of diamond. The parameters for the pressure versus frequency 
relation, Eq. [5] are given. Related experimental data are 
listed in the last two rows of the Table. Note that the lin- 
ear pressure coefficient and the parameter &o are related by 
(du/dP) = too/bo = jo^o/B . 



Method or 


u> 


bo 


V 


(duj/dP) 


Source 


cm -1 


GPa 




cm-VGPa 


All GGA 


1290 


429(6) 


0.045(50) 


3.00 


All LDA 


1322 


460(3) 


0.090(20) 


2.87 


Hanfland et al. a 


1341 


480 


0.118 


2.79 


Nielsen 6 


1306 


456 


0.066 


2.86 


Experiment 


1333 


460 




2.90(5) c 


Experiment 


1333 


446 




3.00(10) d 



"Ref. PL 
6 Ref. [m. 
c R.efs. ^2ql . 

^Corrected value based on revised linear pressure coefficient of the 
ruby Rl line shift. 



of earlier phonon calculations^*^ the P{oj) behavior was 
parametrized using the analog of a Birch expression^ 

P(X)^^b (X 7 -X 5 )[l + r 1 (l-X 2 )} (9) 

with X = (uj/ujo) 1 / 3 . Equation El (more or less an ad 
hoc choice in the earlier work) happens to yield excellent 
representations of the present calculated PJm>) results. 
To facilitate direct comparison with Refs. Q and |27| . 
the parameters &o and r\ are summarized in Table [H] 
Note that for our calculated data r\ is the only adjustable 
parameter, because &o = -Bn/7o- 

Applicability of Eq. 01 in combination with 7 w 1 
implies a nearly quadratic dependence of pressure on 
change in phonon frequency. Thus, we compare calcu- 
lated and experimental results in terms of reduced coor- 
dinates P/Au and Aw, cf. Fig- El Most of the calculated 
results represented in Fig. EI i- e -, the present ones within 
LDA and GGA and those of Refs. Q and j2]j , agree very 
well with respect to the slope. This is of course reflected 
in the small scatter of the r\ values given in Table ITT1 The 
average value from our calculations is 

r) = 0.068 ±0.03 . 

Inserting this value and the experimental data for ojq and 
b Q = (j (dP/duj)o (Table into Eq. yields the P(oj) 
dependence consistent with experimental data near am- 
bient pressure and with the nonlinear behavior predicted 
by the present calculations. 

In the reduced coordinates of Fig. [5] the calculated 
P{ui) results of Wu and Xu^i do not agree well with those 
of other calculations. They use the quadratic function 

lo{P) = lu q + ai P + a 2 P 2 (10) 

to fit their data. The choice of pressure as the indepen- 
dent variable in a quadratic expression is not appropriate 




Frequency shift Aco (cm" 1 ) 

FIG. 5: Calculated pressure as a function of optical phonon 
frequency shift for diamond. A pressure range of about 
200 GPa is covered. Pressure is divided by the frequency 
shift in order to illustrate the nearly quadratic dependence on 
frequency shift. In this representation the intercept at zero 
frequency shift corresponds to the inverse linear pressure coef- 
ficient {duj/dPy 1 . The open circle refers to the corrected ex- 
perimental value. Besides the present LDA and GGA results, 
the figure shows the calculated results of Hanfland et al.r* 
Nielsen)— and Wu & Xu4- The experimental data of Occelli 
et air* are shown for comparison. 



in view of Eq.0|and 7 sa 1. It leads to parameter correla- 
tion which is the possible reason that their results exhibit 
curvature in Fig. 

Occelli et al£ have measured the optical phonon fre- 
quency for pressures up to 140 GPa. They also give 
frequency as a quadratic function of pressure, i.e, a\ = 
2.83 cm-VGPa and a 2 = -3.65 x 10~ 3 cm^/CPa 2 in 
Eq.EH The curvature of their data in Fig. \5\ could again 
be related to the particular choice of the fitting func- 
tion. The main observation, however, is that the aver- 
age slope for the experimental data differs from the one 
predicted by theory. At the frequency corresponding to 
140 GPa experimental pressure the pressure difference 
between calculations and experiment again amounts to 
about 11 GPa. 

The sign and magnitude of this difference is very sim- 
ilar to that encountered when comparing experimental 
and calculated EOS results. Thus, in terms of pressure 
we find differences between calculated and experimental 
results which are about the same for different physical 
quantities considered. This hints to an explanation where 
the discrepancies are caused by the same systematic er- 
ror (s). 
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V. REMARKS ON THE RUBY CALIBRATION 

We consider the possibility that the discrepancy in the 
high pressure regime between experimental and calcu- 
lated EOS and phonon frequency results for diamond are 
caused by some error in the experimental pressure scale, 
i.e. the calibration of the ruby Rl line shift according to 
Ref. 

Holzapfe! 13 recently proposed a revised ruby pressure 
scale, based on an analysis of published EOS data of se- 
lected elemental solids, including those of Ref. yj for dia- 
mond. The revised ruby calibration was cast into a three- 
parameter analytical expression for pressure as a function 
of the Rl line wavelength. With the recommended pa- 
rameter values^ and when restricted to about 200 GPa 
pressure (a range beyond the upper limit of most DAC 
studies), the revised calibration can be represented by a 
simple second-order polynomial in frequency shift Av of 
the Rl line. This is evident from Fig.[S]which shows plots 
of (| A^l/t'o) 1 x P versus \Av\/vq for Holzapfel's revised 
scale and for the 1978/1986 scales of Mao et alM& We 
write the quadratic polynomial asi& 



P{Av) = A 



B |Ai 



2 v Q 



(11) 



A linear regression for the 'Revised' line in Fig. [S] with 
A = 1820 GPa (fixed value as recommended in Ref. 01) 
gives a slope parameter of i?=15.8. 

From the experimental P(V) and P(oj) data of 
diamondZ one can recover the ruby line shift Av accord- 
ing to the ruby calibration^* used in the experimental 
work. Combining the experimental V(Av) and u(Ai>), 
respectively, with the calculated P(V) and P(co) rela- 
tions reported here yields P(Av) as shown by the two 
solid lines in Fig. El These curves run more or less par- 
allel to Holzapfel's revised calibration. Obviously, both 
the calculated P(V) and P(u>) relations for diamond in 
combination with the experimental data of Ref. Q sup- 
port the proposed revision of the ruby pressure scale. It 
should be noted that the ruby calibration discussed by 
Aleksandrov et al. 5 would be a little higher in pressure, 
but in the coordinates of Fig. it exhibits a slope similar 
to that of the 'revised' line. 



VI. CONCLUSIONS 

We offer the following conclusions: 

(1) The results of the first-principles EOS and phonon 
frequency calculations for diamond reported here do not 
depend on the computational method, i.e., the choice of 
the crystal potential and basis sets in the PW, PAW, and 
APW+lo methods. 

(2) The calculated equilibrium values for volume (Vq) 
and bulk modulus (-Bo) of diamond do depend on the 
exchange-correlation functional (LDA or GGA), a well- 
known fact in density functional theory. However, the 



3000 



2800 



7a 2600 

Q_ 

(5 

J 3 " 2400- 



3 



2200- 



2000 



1800 



oc-AI 2 3 :Cr 

R r line 



Revised 



Diamond phonon / 
/ 

/ 

/ 

/ 

/ / 1 986 

Diamond EOS 




''1978 



MgO elast. 
Piermarini, Nakano 



- Grasset R-,, R 2 

I L 



0.00 0.02 0.04 0.06 0.08 
Energy change |Av|/v 

FIG. 6: Ruby pressure calibrations in reduced coordinates. 
Pressure is divided by energy change of the Ri line and plot- 
ted as a function of the energy change. Dashed lines refer to 
the calibrations by Mao et alw—Z— (marked 1976 and 1986) 
and the revision proposed by Holzapfel^ Lines marked 'dia- 
mond' refer to converted data of Occelli et alw- combined with 
the calculated EOS and phonon frequency shift of diamond 
reported in the present work. The symbols near zero energy 
shift correspond to Rl line pressure coefficients determined 
at low pressures!— The symbol marked 'MgO' stands for 
the result of elasticity studies of MgO up to 55 GPai— 



value of the pressure derivative of the bulk modulus 
[-B =3.65(5)] and the analytical form of the P(V) relation 
are not affected within the volume range (40% compres- 
sion) covered here. In other words, the calculated non- 
linear component in the P(V) behavior is independent of 
the method of calculation and the exchange-correlation 
approximation. 

(3) The different theoretical methods and approxima- 
tions also yield very similar results for the nonlinearity 
in the pressure versus phonon frequency relation P{lo). 

(4) Both the calculated P(V) and P(u>) relations ex- 
hibit a nonlinear behavior which is more pronounced than 
that observed in recent x-ray diffraction and Raman mea- 
surements of diamond up to 140 GPai In terms of ab- 
solute pressure, the differences between theoretical and 
experimental results are essentially the same for the P{V) 
and P(oj) behavior, of the order of 10 GPa at an experi- 
mental pressure of 140 GPa. 

(5) Our calculated nonlinearity of P(tS) agrees well 
with that obtained in earlier ah initio phonon frequency 
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calculations for diamond under pressureASi The work of 
Hanfland et alA indicated the need to reconsider the cal- 
ibration of the ruby pressure scale used at that timeiS 
The 'quasi-hydrostatic' (1986) ruby calibration 14 reduced 
the discrepancy between theoretical predictions and ex- 
perimental results. The remaining difference discussed 
here at least calls for some caution when experimental 
data measured above 50 GPa and based on the 1986 ruby 
pressure scale are compared to ab initio calculations or 
reduced shock wave data. 
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